Topological confinement in bilayer graphene 
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We study a new type of one-dimensional chiral states that can be created in bilayer graphene 
(BLG) by electrostatic lateral confinement. These states appear on the domain walls separating 
insulating regions experiencing the opposite gating polarity. While the states are similar to conven- 
tional solitonic zero-modes, their properties are defined by the unusual chiral BLG quasiparticles, 
from which they derive. The number of zero-mode branches is fixed by the topological vacuum 
charge of the insulating BLG state. We discuss how these chiral states can manifest experimentally, 
and emphasize their relevance for valley tronics. 



Most condensed matters systems acquire gap in single 
electron excitation spectrum at low temperatures. Typ- 
ically this happens as a result of spontaneous symmetry 
breaking, as in the case of superconductors or charge- and 
spin-density-wavc materials. The gap opens due to the 
interaction of electrons with a slow bosonic mode repre- 
senting the order parameter (OP). Yet, for certain topo- 
logically non-trivial configurations of OP, zero-energy 
fcrmionic states reemerge. Examples are Andreev states 
that form at the domain walls in superconductors [H, 
states in the superconducting vortex cores 0], and soli- 
tons in polyacethylene Q . Similarly, in cosmology, it has 
been suggested that our 3+1 dimensional space with its 
extremely low elementary particle masses may represent 
a membrane or a string in a higher-dimensional inhomo- 
geneous Higgs vacuum The zero modes may exhibit 
a number of interesting phenomena, including fermion 
number fractionalization [5| and chiral anomaly [6J. In 
condensed matter, similarly to cosmology, the zero modes 
originate from the Dirac equation, which emerges as an 
effective (linearized) description of physics near the Fermi 
surface. 

In this Letter we demonstrate that zero-modes can also 
emerge in electrostatically-gated bilayer graphene (BLG) 
structures (Fig. [1]). These modes, however, are different 
from all the examples from condensed matter and high- 
energy physics that we know of, as they derive not from 
the Dirac fermions [5|, M] , but from the unusual low-energy 
chiral modes of BLG [8j , which have quadratic dispersion 
and zero gap between particle and hole bands. When 
electrostatic bias V is applied between the layers, a gap 
of the magnitude \ V\ opens between the particle and hole 
bands 0. The interlayer bias plays a role analogous to 
an OP, but is externally tunable. By applying inhomoge- 
neous bias one can spatially confine the low-energy states 
to the regions with low gap, e.g. to a one-dimensional 
(ID) channel. However, in addition to the conventional 
confinement by constant-polarity potential, in the case 
of BLG, there is a possibility of topological confinement, 
with the sign of the confining potential changing across 
the channel, Fig. [T] We find that the topological con- 
finement leads to formation of ID chiral zero-modes. In 



each valley of graphene band structure there are two such 
modes (per spin) , both carrying electrons in the same di- 
rection (opposite for the two valleys). The robustness 
of the zero-modes and their chiral nature are ensured 
by the topological structure of the gapped bulk states. 
These modes are likely to influence transport in BLG 
with smooth potential disord er [ldl| . They may also have 
implication for valleytronics [111 ], as they can enable the 
fabrication of valley filters and valves, which can be ex- 
perimentally realized with existing technology. 
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FIG. 1: Side view of gated bilayer graphene configuration with 
the voltage kink. The region where the interlayer voltage 
changes sign (channel) supports bands of chiral zero-modes 
(dashed line). The conventional (non-topological) confine- 
ment would correspond to the same polarity of bias on both 
sides of the channel. 



The low-energy (one-valley) bilayer Hamiltonian with 
bias V(x) applied between the layers is d, @] 
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where c is the Fermi velocity, it = p x + ip y , tt^ = 
Px — ipy, and t± is the interlayer coupling. The Hamilto- 
nian acts in the space of four-component wavefunctions, 
(ipAi)' l ^Bi," l pA2) ^52), where letters in subscript represent 
sublattice and numbers - the layer (we included here only 
one valley, assuming that V(x) is smooth enough not to 
cause intervalley transitions). This Hamiltonian provides 
a good description of BLG as proven by recent experi- 
ments 12 1. At low energies and constant gate voltages, 
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V <C t± , the Hamiltonian can be further reduced [lc 
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Here p = y P 2 + P y ■ The remaining components of the 
wavefunction have predominantly Al and B2 character. 
The spectrum of the Hamiltonian is 
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At finite bias, the spectrum has a gap which reaches min- 
imum |V| on the circle cp ~ V/2. However, as long as 
V < t±, the p 2 correction in the first term can be ne- 
glected, simplifying spectrum to E 2 — V 2 /A + c p /t\. 
The corresponding term can be also neglected in the di- 
agonal part of H. As a result, in the Hamiltonian (0), 
the V terms and the momentum terms are decoupled. 
Thus, the transition to the position-dependent potential 
V(x) corresponds to reinstating the momenta p x and p y 
as differential operators, finally leading to a dimension- 
less quasiclassical Hamiltonian [14j |. 

T-tqc = -f{x)a z -{p 2 x -p 2 y )a x -2p x p v ay = g(p,a;><7, (4) 

where we defined <p{x) — V(x)t±a 2 /(2c 2 ), and momenta 
are measured in units of inverse lattice constant a. 

For constant ip(x) the spectrum has the form E = 
±-\/A 2 + e(p) 2 , which is similar to the spectra of con- 
densed matter systems with off-diagonal long range or- 
der; the gate voltage ip plays the role of the order param- 
eter A. It is well known that in such systems with non- 
trivial topological structure of A(r), low-energy fermionic 
modes can emerge. Motivated by this analogy, we study 
related inhomogeneous configurations of ip{x). 

From Eq. (j4]), the corresponding wave equation is 



(p(x)u + (d x +p y ) 2 v = su, 
tp(x)v + (d x - p y ) 2 u = ev. 



(5) 
(6) 



It possesses a number of symmetries for a general an- 
tisymmetric potential profile, ip(—x) — —ip(x). It is 
easy to see that in this case for v(x) = ±u(—x) the 
two equations (0) and (0) reduce to one. Therefore one 
can solve two problems separately, for *S? = [u(x),u(—x)] 
and $ = [w(x), —w(—x)]. Furthermore, after we find the 
eigenvalues and eigenvectors of the first problem for some 
value of p y , = [u™ v (x),u™ y (-x)} and e£ v , the other 
solution can be obtained as <&™ y = [— u" Ph (— x), u" Ph (x)) 
with eigenvalues (— e" p ). Thus in the p y — e plane the 
dispersion of $ is obtained from dispersion of "J" by inver- 
sion relative to the point (0, 0). A related, useful symme- 
try has to do with the behavior of the spectrum of Eqs. 

and ^ under p y — > — p y . Easy to see that under this 
transformation (u, v) — > (v, —u) and e — > —e. 



Unlike the ID Dirac equation, which typically de- 
scribes the domain- wall zero-modes [5], Eqs. J5]) and (J6]) 
are not easily solvable analytically for a general profile 
(f(x). We therefore first analytically study the simple 
case of step-like potential ip{x) = ipo sign(a;), and then 
perform numerical solution of Eqs. ([5]) and © for a 
smooth ip{x). As we will see, only the details of the 
zero-mode dispersion e(p y ) depend on the exact profile 
of ip(x). This is a manifestation of the topological nature 
of these states. 

Step kink. We first consider a step-like kink, <f(x) — 
</7osign(x). In this case, both for positive and negative 
x the potential is constant, and the solution of the wave 
equation in these regions is 9 oc e~ Xx . For the intragap 
states, i.e. those with |e| < ipo, 



A = ±q ± i(3, 



(7) 



where a(fi) = 2- l / 2 [{p y + V 2 -e 2 ) 1 ' 2 + (-)p 2 y } 1/2 . For x < 
one should only keep A^ 2 = a ± ij3 and for x > only 
\f 2 = —ocH/3, with the corresponding wavefunctions of 
the form 

u^(x) — uf exp(— Xf x) + uj- exp(— AJx). (8) 

At x = the left and right solutions have to be matched. 
From the structure of the wave equation it is clear that 
the wavefunction and its first derivative are continuous 
across x = 0, while the second and higher derivatives are 
not. Considering for concreteness the states of the form 
= [u(x), u(— x)], the matching conditions are 
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The third equation is obtained by subtracting Eq. ([5|) at 
x = — from itself at x — +0 and using v(0) — u(0). 
The forth equation is obtained in the same manner but 
after first differentiating Eq. ([5]) over x and using that 
d x v(0) = —d x u(0). Substituting now the general solu- 
tions Eqs. ijHJ), we obtain a homogeneous system of 4 
equations with 4 unknowns which only has a non-trivial 
solution if its determinant equals zero. Given the form 
of A, Eq. |0), this condition is equivalent to 
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Near zero energy, this equation has a solution only for 
p y < 0. For (fo > 0, we obtain the dispersion 
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(we analyzed the wave equation near e = to remove 
a spurious branch with e — > — e). The zero-energy solu- 
tion obtains at po = — y^o/2 3 ' 4 . The other branch cor- 
responding to wavefunction $ = [— u- Py (— x), u_ p (x)] 
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has, as discussed above, the inverted dispersion, 

Pv = (-e + ip /V2)/(-e + <p V2) 1 / 2 . (15) 

Notice that both solutions have negative velocity near 
e = 0. This seems to imply time-reversal symmetry 
breaking. However, for the second graphene valley, the 
velocity is positive and thus the symmetry is reinstated. 
Therefore, the zero-modes are chiral: i.e. if we define 
pseudospin-1/2, with S y corresponding to the valley in- 
dex, all zero-modes have definite sign of p y S y . It is also 
easy to see that on the anti-kink (</jq < 0) the dispersion 
is flipped, e{p y ) — > -e{p y ). 

Note that the decay length of the wave function 
of the topologically confined states is of the order of 
a/y/ipo pa at/^/Vt± > a at low energies. Therefore, our 
solution is consistent with the long-wavelength expan- 
sion (|4]) used for the description of the system. For kinks 
wider than a/^/Ipo one expects quantitative deviations, 
which we indeed find in the direct numerical solution. 
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FIG. 2: Spatial structure of the confined state obtained nu- 
merically from Eqs. ((5| and <(6j for particular value p y — 
— 1/2 3/4 . The energy of the state is e = 0.1140. Shown are 
the components it and v of the wavefunction (thin and dashed 
lines, respectively), the probability density |^| 2 = u 2 + v 2 
(thick line), and the potential profile, taken here as <p(x) = 
tanh(x). Note the symmetry, u(x) = v(— x). 

Numerical solution. The wave equations ([5]) and ^ 
can be solved for an arbitrary potential profile if(x) nu- 
merically using the finite differences method. In Figure 
[2] we show an example of an electronic state localized at 
the kink of the potential ip = tax\h(x/£) for i = 1 and 
p y = — 1/2 3 / 4 . For step-like potential this value of p y 
would correspond to a zero-energy state. For the rather 
smooth potential used here, the lowest energy is finite, 
e = 0.1140, but indeed much smaller than the reference 
gap, tpo = 1. The state has the expected symmetry, 
v(x) — u(—x). The wavefunction has a slightly oscilla- 
tory and decaying behavior, consistent with the complex 
(non-real) values of A's. A symmetry-related low energy 
state $ = (— v,u) occurs at p y = + 1/2 3 / 4 . 




FIG. 3: Energy level structure in the presence of non-uniform 
interlayer bias ip(x). The low-energy states in (a) and (c) are 
localized on the single kink of profile <fi{x) = tanh(a;/i?) with 
I = 0.5 and £ — 4, respectively. For (a) and (c) we used open 
boundary conditions, ^(±20) = 0, which caused flat bands at 
e = ±1. The dashed lines are the analytical expressions for 
the intragap state dispersion, Eq. 1(140. (b): Same as (a) but 
with kink and anti-kink at x = ±10 and periodic boundary 
conditions, (d): Non-topological confinement for potential 
profile ip(x) = tanh 2 (a;/4). Note the absence of chiral zero- 
modes. 



In Figures^-c we show the intra-gap state dispersions 
e(p y ) for kinks of various widths. Note that for the nar- 
row kink (panels a and b) there are only two intragap 
states, while for the wide kink (panel c) more low-energy 
states appear below the gap edge. In all cases, there are 
only two (chiral) zero-modes per kink per valley per spin 
regardless of the width. As expected, the analytical re- 
sult, Eqs. (fl4|) and (fl~5|) , fits very well for the narrow kink. 
For comparison, panel d shows the dispersion relation for 
states in a ID channel defined by "conventional" confine- 
ment. In this case no chiral zero modes are present, even 
though there are states below the gap edge. 

Charge of the kink. It is well know that solitons in ID 
systems can carry char ge, which can be either rational 
or irrational 

@, El EUl El- Similarly, m our case, 
one might expect that the kink can carry charge. Let 
us demonstrate that contrary to this expectation the do- 
main wall charge density is zero. Using the symmetries 
of the problem, the density in the presence of a kink is 
proportional to 
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By completeness property of the wavefunction, this ex- 
pression is an ^-independent constant proportional to the 
momentum (energy) cut-off. 

Topological considerations. Above we found that in the 
presence of a kink there are two zero-modes per graphene 
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valley per spin. We now demonstrate how the appear- 
ance of these modes can be understood from topologi- 
cal considerations (l9| . The number of zero- modes, or 
more precisely the number of zero-energy left-movers mi- 
nus the number of zero-energy right movers is related 
to the topological charge A3 of the Fermi point located 
in the extended quasiclassical 3D momentum-real space 
(p,x). In terms of the quasi-classical Green function, 

Gqc = (iPO - H«j C )~\ 
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where integration is to be performed over an arbitrary 3D 
surface 03 enclosing the Fermi point (po, p, x) — (0, 0, 0). 
The charge Na is the degree of the mapping of 0-3 onto 
a manifold corresponding to the Q qc . Given the form 
of Tl qc , this manifold is equivalent to a 3D sphere, and 
thus the mapping belongs to non-trivial ira(S 3 ) homo- 
topy group. In evaluating Na, two choices of 03 are par- 
ticularly useful: (a) infinite planes p y — ±p® and (b) 
infinite planes at x = ±x° (assuming that the planes 
are sufficiently close, the flux contributing to Na but 
not crossing these planes is negligible). The represen- 
tation (a) is equivalent to the quasiclassical expansion 
of the difference of spectral asymmetry functions, v — 
v (P°y) - v{~Py), wh ere 
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v counts the number of dispersion branches that cross 
zero energy from above on the interval [~Py,Py]- In the 
representation (b), A^ is nothing but the difference of 
the vacuum topological charges of the insulating states 
to the right and to the left of the domain wall, A^ir ) — 
N 3 (-x°). Given the form of H qc , Eq. |@J), 

If 1 

N s(xo) = I dp x dp y j^g ■ [dp^g x d Py g], 



from which one easily finds Na(xo) = sign[<yj(:r )]. Thus, 
for <p(+oo) > 0, we obtain N 3 = N 3 (x ) — N 3 (—x ) = 
v = 2, consistent with our earlier finding that there are 
two branches per kink per valley per spin that cross zero 
energy with negative velocity. 

Having established the presence of one-dimensional 
chiral zero-modes originating from topological confine- 
ment, we briefly discuss their relevance to one area of 
current interest, namely " valley tronics," which attempts 
to utilize the valley de gree of freedom to achieve new 
electronic functionality [llj. As we have shown above 
(e.g. see Fig. [3£i,c), a topologically confined ID channel 
contains zero modes, with the direction of propagation 
determined by the valley. Thus, a valley filter is realized 
when a voltage difference is applied along the kink - only 
electrons in one of the valleys will carry the current and 



the other valley is "filtered out." It is also straightforward 
to realize a valley valve by simply "connecting" two filters 
in series. When the polarity of two filters is the same, 
the current can pass through; for opposite polarities the 



current is blocked. We note that earlier proposals [111.12 



for valley filters and valves in graphene monolayers relied 
on perfect zig-zag edges, making their practical realiza- 
tion very challenging. On the other hand, the technology 
needed for valley filter and valve using the topologically- 
confincd channels is much less demanding: A smooth ID 
channel can easily be realized away from the edges of the 
graphene sample, thus decreasing substantially the pos- 
sibility of intervalley scattering. Some of the technology- 
the opening of a gap in a bilayer using gate electrodes-has 
already been demonstrated experimentally 2l|. A more 
detailed discussion of valley filtering induced by topolog- 
ical confinement will be presented elsewhere 22j. 
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